Semiconductor quantum tubes: dielectric modulation and excitonic response 
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We study theoretically the optical properties of quantum tubes, one-dimensional semiconductor 
nanostructures where electrons and holes are confined to a cylindrical shell. In these structures, 
which bridge between 2D and ID systems, the electron-hole interaction may be modulated by a 
dielectric substance outside the quantum tube and possibly inside its core. We use the exact Green's 
function for the appropriate dielectric configuration and exact diagonalization of the electron-hole 
interaction within an effective meiss description to predict the evolution of the exciton binding energy 
and oscillator strength. Contrary to the homogeneous case, in dielectrically modulated tubes the 
exciton binding is a function of the tube diameter and can be tuned to a large extent by structure 
design and proper choice of the dielectric media. 

PACS numbers: 73.22.Lp, 78.67.Ch 
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I. INTRODUCTION 

Cylindrical semiconductor nanostructures bridge be- 
tween quasi-lD systems at small diameters and quasi-2D 
in the opposite limit, thus extending the wealth of physics 
and applications of low-dimensional solid-state systems. 
The controlled growth of semiconductor quantum tubes 
(QTs) with diameters in the 10-100 nm range has been re- 
cently demonstrated through several techniques, includ- 
ing multi-layer overgrowth of nanowiresi"— and strain- 
induced bending of a planar heterostructurei^'^ In addi- 
tion to QTs with a solid semiconductor core, it is possible 
to grow hollow QTs, where the charge carriers are con- 
fined in a thin semiconductor shell, encompassed by a 
barrier material which is only a few nm thick.i^"^ Large 
surface-to-volume ratios and the possibility of various 
functionalizations on both the internal and external sur- 
faces make the latter systems particularly interesting for 
applications!^ 

Although experiments concerning the optical proper- 
ties of these systems are still limited, advancements in 
the optical quality of the samples point to a rapid in- 
crease of these investigationsi^iiSrJ^ The excitonic prop- 
erties of semiconductor QTs are particularly interest- 
ing with respect to conventional semiconductor quan- 
tum wires, where excitons are confined in the core of the 
nanostructure.— On the one hand, due to the com- 
bined effect of the QT curvature and of the quasi- 2D 
confinement of carriers in the cylindrical shell, excitonic 
binding energies might be substantially stronger than in 
bulk, even for large diameter QTs. On the other hand, 
a dielectric medium outside the shell of the QT may re- 
sult in a dielectric confinement of the electric field felt 
by the optically excited electron-hole pairs, in most cases 
enhancing their excitonic binding energy. Since the di- 
electric interface is spatially separated from the carriers, 
which are confined deep inside the shell, excitonic bind- 
ing and sensitivity to the medium might be strongly en- 
hanced without spoiling the optical properties of the elec- 
tronic system;^ analogously to core-shell nanowiresJ^ 



The screening provided by the dielectric environment can 
be varied in a broad range<2^ The tunability of the di- 
electric constant in the core of the QT, obtained, e.g., by 
oxidation;^ can further increase such effects. 

Present work on QTs theoretically considered mag- 
netic statesr^i^ reported experimental evidence of 
the Aharonov-Bohm effect,— and treated optical 
properties;^ but the influence of the dielectric dismatch 
between the nanostructure and the environment has been 
studied so far only for conventional quantum wires^i and 
freestanding nanowiresi^ Here we will consider also a di- 
electric mismatch between the core and the shell, which 
will lead to a considerable change in the electron-hole 
interaction, as shown in Fig. [2l 

Hereafter we investigate the excitonic binding and os- 
cillator strength in hollow and filled QTs for different 
geometries and dielectric configurations. Besides increas- 
ing due to the reduced screening, the excitonic binding 
strongly depends on the QT diameter and on the dielec- 
tric medium. The paper is organized as follows. In Sec- 
tion In] we outline the theoretical model, which includes 
the exact solution of the Poisson equation and the diag- 
onalization of the electron-hole Hamiltonian within the 
envelope-function approximation. In Sections HI and IV 
we report our results and draw the conclusions, respec- 
tively. 



II. THE MODEL 

The system we consider consists of an infinite tube with 
cylindrical symmetry^^ (see Fig. [1}. For simplicity we 
assume that the motion in the radial direction is frozen, 
and that charge carriers are radially confined in a 5-like 
well at a distance R from the tube axis. This electronic 
layer is buried in the middle of a coaxial cylindrical shell 
of thickness 26 with dielectric constant £5, while the core 
and the environmet have in general different dielectric 
constants, ec and es, respectively. Since the shell is a 
semiconductor material, typically eg > ec,e-E^ 
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FIG. 1: (color online) Schematics of a QT with outer surface 
cut open. The electron-hole pair are constrained to a cylin- 
drical surface of radius -R, embedded in a shell of thickness 
25. X and y are the relative electron-hole coordinates. The 
dielectric constants eE,€s,^c characterize the environment, 
the shell and the core region, respectively. 



The invariance under translations along, and rotations 
around the tube axis warrants the separation of the cen- 
ter of mass and relative coordinates. The motion of the 
Wannier exciton^ in the relative degrees of freedom is 
determined by the envelope-function Hamiltonian 



dx 



(1) 



expressed in units of the effective Hartree Ha* = 
(/i/e|)Ha, with ^ = merrih/ (rrie + ruh) the reduced 
electron-hole mass. The relative coordinates around the 
circumference {x — R(j)) and along the tube axis (y, 
see Fig. [1]), are in units of the effective Bohr radius, 
= (es/M)aB. 

The effective Coulomb interaction potential V{x, y) be- 
tween the confined electron and hole depends paramet- 
rically on the dielectric constants (ec eg, e^) and on the 
tube geometry through i5 and R. In cylindrical coordi- 
nates, the potential (scaled with Ha*) generated by a 
charge at r' = {p' , (j)' , z') reads 



dkcos[k{z - z')\gm.a{k,p,p'), (2) 



where a = C, S, E indicates whether the position of the 
test-charge r = (p, </>, z) is in the core, shell or environ- 
ment region, respectively, and gm,a{k, p, p') is the solu- 
tion of the radial Poisson equation in that region (see 
the Appendix for further details). The interaction V in 
Eq. ([T|) coincides with Vs, with p = p' = R. As shown in 
the Appendix, 
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gm,s{k,R,R) = — 



ikR)K„,{kR), 



CI 



where Im , Km are the Bessel functions of the first and 
second kind; the coefficients B^,C^,B^,C^ are given 
in Eqs. (|A.10[) and (jA.12[) in terms of Im, Km and their 
derivatives. 

To illustrate how the electron-hole interaction is infiu- 
enced by the dielectric environment, we shown in Fig. [2] 
the potential V for i) a filled QT, with a core of the 
same material as the shell, immersed in a substance 
with a low-dielectric constant (ep = es = lOe^), and 
ii) a hollow QT, with the same low-dielectric constant 
substance inside and outside the shell (eg = lOe^; = 
lOes. For comparison, we also show the dielectri- 
cally homogeneous case (ec = es = ^e), where the 
V reduces to the usual Coulomb potential V{x, y) = 
-l/esV(2-Rsin(a;/2i?))2 -|-y2^ FigurelSJa) shows the in- 
teraction along the QT V{x = 0, y), while Fig.[2jb) shows 
the interaction around the cylinder V{x,y = 0). The 
Coulomb interaction for the hollow and filled cases is for 
all distances stronger than in the homogeneous case, since 
the average dielectric constant of the system is smaller, 
and the electric field is not screened outside and, for the 
hollow case, also inside the QT. The interaction in the 
filled and hollow case is substantially different only for 
distances smaller or comparable to the Bohr radius, with 
the interaction in the hollow case being stronger. For 
larger distances (inset of Fig. H^a)), on the other hand, 
the non-trivial influence of the dielectric mismatch be- 
tween the core and the shell leads to crossing of the po- 
tentials for hollow and filled QTs, before both converge 
to the same value. 
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FIG. 2: Electrostatic interaction V{x,y) between an electron- 
hole pair confined to a cylindrical surface with diameter D = 
0.8 a^, buried in a shell with thickness 25 = 0.1 a^ with 
dielectric mismatch, as follows. Homogeneous case: ec = 
es = (-E. Filled case: ec = ts = lOtB. Hollow case: es = 
lOec = lOtB. (a) interaction along the QT, V{x = 0,y). 
(b) interaction around the QT, V{x,y — 0). Inset in (a): 
V{x = 0, y) in a larger range of y. 



A convenient basis set to represent the cxciton wave- 
function is obtained by multiplying eigenfunctions of 
the linear momentum operator along y (e**"'^) and of 
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the angular momentum operator along the tube axis 
(e™^/^). Imposing periodic Born- von Karman bound- 
ary conditional along y, with period L sufficiently larger 
than the effective Bohr radius of the material, results in: 
k = pAk with Afc = 2t: / L. The wavefunction thus reads 

N P 

E E C';i^pe--/«e'f^^■^ (4) 

n=-N p=-P 

where p,n Cz and j indicates the j— th exciton state. 
The coefficients ^ are obtained from the Schrodinger 
equation in the above basis: 



u„ 



r<3 — wci 



(5) 



The diagonal term in the first line represents the kinetic 
energy, whereas the matrix elements of the electron-hole 
Coulomb interaction term are given by 



Afc 



g\n-n'\,s{\p-p'\Ak,R,R). (6) 



In order to reduce the dimension of the Hamiltonian ma- 
trix we introduce a cutoff energy Ecut, set the maxi- 
mum number of plane waves P = V2£'cut/ Afc and choose 
the maximum number of orbital modes N in Eq. Q as 
the nearest integer to n{p) = R^y2Ecut — {pAk)^. The 
Hamiltonian matrix is block diagonalized using a sym- 
metrized basis set. In particular, we consider linear com- 
binations of the above basis functions that are even or 
odd with respect to the inversion of the relative coordi- 
nates X and y, which is the equivalent of inverting the ab- 
solute coordinates, since the corresponding inversion op- 
erators Tlx , Tly commute with the relative motion Hamil- 
tonian 







0. 



(7) 



The resulting energy Ej is obtained with respect to the 
energy minimum of the conduction band. Therefore the 
binding energy of the exciton ground state is E^ = —Eq. 
In the presence of a photon gauge field the electron-hole 
pair recombines emitting a photon of energy Eg — Ef. 
The recombination rate is related to the dimensionless 
oscillator strength /, which in the dipole approximation 
reads^ 



k'o(0,0)|% 



E„ 



(8) 



Here |'0o(O,O)p is the envelope function of the exciton 
ground state given by Eq. (|4]), Q the momentum of the 
center of mass. Eg is the energy gap between valence and 
conduction band and 5o — Ep/tg, where Ep is the energy 
associated with Kane's matrix elementsi^ 



III. RESULTS 

In the following we investigate the excitonic proper- 
ties of QTs made of the direct gap materials, InAs, 
GaAs and InP, and two different dielectric configurations: 
filled QTs, with a core of the same material of the shell 
{ec = es ^e), and hollow QTs, with the core of the 
same material as the environment (ec = 7^ £5). We 
consider QTs with diameters in the 20 100 nm range 
and a constant shell thickness 2S = 10 nm, compara- 
ble to state-of-the-art samplesjiiii^ Material parameters 
used in the calculations are listed in Table HI 

In Fig. [3] we plot the energy of the exciton for hollow 
and filled QTs in vacuum (eg = 1), and compare it with 
the 2D limit i? — > 00; this, in the case of excitons confined 
to a strictly 2D layer, is four times the bulk valuei^"— 
The exciton binding energy shows a strong increase with 
respect to the 2D limit, and a marked diameter depen- 
dence, which is different for hollow and filled QTs: while 
in the former case the exciton binding energy becomes 
weakly dependent of the diameter for QTs larger than 
~ 60 nm, the latter one show a strong dependence even 
for the larger QTs. 

It is instructive to contrast these results with the ex- 
citon energies of QTs which are dielectrically homoge- 
neous, that is, buried in a material with the same bulk 
dielectric constant of the semiconductor shell (es = ec = 
e^;). In this case, exciton energies do not show any depen- 
dence on the diameter, and are pinned to the 2D value 
(see Fig. [3]) . This is due to the small value of the Bohr 
radius with respect to the tube diameters shown here, 
so that the curvature of the surface plays only a minor 
role.— Clearly, for smaller diameters (not shown here), 
the exciton binding energy increases and the exciton en- 
ergy red-shifts, since the binding energy is infinite in the 
strictly ID limitj^i^ In Fig.[3]this can only be recognized 
in the tiny red-shift for InAs at the smallest diameter. 

The large binding energy of dielectrically modulated 
QTs with respect to homogeneous ones is an obvious 
consequence of the smaller screening of the electron-hole 
interaction in regions where the dielectric constant is 1. 
Accordingly, exciton energies of filled and hollow QTs 
are similar for the smaller diameters, because their di- 
electric configuration differs only in the core, which is a 
small fraction of the volume if 26 = 10 nm. Increasing 
the diameter while leaving the shell thickness constant 
corresponds to increasing the core region with respect to 
the shell, thus enhancing the difference between filled and 
hollow QTs. In both cases energies are increasing with 
diameter due to the larger area occupied by the cross- 
section of semiconductor shell, but faster for filled than 
for hollow QTs. This is consistent with the following ar- 
gument: the cross section of the tube with es is a ring of 
area 2D5Tr for hollow QTs and a circle of area {D/2+S)^Tr 
for filled ones. Thus, the screening area is growing faster 
for filled QTs. In fact, the exciton energy for hollow and 



filled QTs increases as 
tively. 



-lO/D and lO/D^, respec- 
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FIG. 3: Ground state energy of the exciton for InAs, GaAs and InP tubes. Different symbols correspond to a homogeneous 
dielectric constant (ec = es = ^e, crosses), a filled tube {ec = es, = 1, filled circles) and a hollow tube (ec = = 1, empty 
circles). The shaded region marks the analytical 2D limit for the homogenous case. The exciton energy for hollow and filled 
QTs increases as ~ — 10/D and ~ — lO/D^, respectively. 



All in all, excitons in QTs can be from twice (filled 
GaAs QTs, D = 100 nm) up to 7 times (hollow InAs QTs 
at Z? = 20 nm) more strongly bound with respect to the 
respective 2D bulk exciton. Moreover, in all investigated 
cases they have a binding energy which is much larger 
than thermal energy at room temperature. 

In Fig. [3] we used the permittivity of vacuum, which 
leads to the largest possible dielectric confinement effects 
in a given configuration. Next we show how the binding 
energy depends on the dielectric constant of the medium 
by which the QTs are surrounded i^i^^i^ Figure S] shows 
the binding energy as a function of the ratio es/^s- The 
leftmost point on the horizontal axis corresponds to the 
case = 1, while the ratio eb/es = 1 correspond to the 
homogeneous case, with the same dielectric constant es 
filling all space. The general behavior is the same for the 
three materials. The difference in energies between hol- 
low and filled tubes is largest in vacuum for large diame- 
ters. Increasing the screening of the surrounding leads to 
less strongly bound excitons, since the Coulomb interac- 
tion is more and more inhibited, and the two cases of hol- 
low and filled tubes are becoming increasingly similar to 
each other, and obviously coincide at e^j/eg = 1. We also 
note that the exciton binding energy is very sensitive to 
the environment (and core) in the low-dielectric-constant 
range. The exciton binding energy with respect to the 
2D case is halved increasing ee from O.les to ~ 0.2e5, 
which suggests that small changes of the dielectric envi- 
ronment might be revealed by optical means in this type 
of system, which is due to the proximity of the electronic 
system to the environment. We stress that excitons in 
InAs tubes are always less bound than excitons in GaAs 
or InP tubes for any considered value of the ratio ee/es- 

In Fig. m we also show in color-scale (color online) and 
point-size the oscillator strength of Eq. ([8]) as a func- 
tion of diameter and dielectric configuration. For each 
material the oscillator strength is normalized relative to 
the maximum value for the same material, so that the 



TABLE I: Parameters of materials under consideration. The 
mass, energy and length are in units of bare electron mass, eV 
and nm respectively. Values are taken from Refs. [3^l4lll42 . 



Material 






es 


Ha* 


a*B 


Eg 


InAs 


0.026 


0.33 


14.6 


0.0031 


32.05 


0.35 


GaAs 


0.067 


0.35 


12.5 


0.0098 


11.76 


1.43 


InP 


0.08 


0.33 


12.4 


0.0123 


9.23 


1.34 



atomic part of the oscillator strength cancels out. The 
results shown in Fig. |4] can be summarized as follows: i) 
For GaAs and InP QTs the recombination probability is 
larger the smaller is the diameter, while for InAs QTs it 
is nearly insensitive to it; ii) for GaAs and InP QTs the 
relative oscillator strength is very weakly dependent on 
EE, while for InAs QTs it decreases for increasing ratio 
^e/^s- The peculiar behaviour of InAs with respect to 
GaAs and InP in both respects must be traced to inter- 
play between the small band gap of InAs and the very 
large exciton binding energy in this class of systems, mak- 
ing the denominator in Eq. ([5]) strongly dependent on . 

In order to further investigate the effect of the Coulomb 
interaction between the carriers, we plot the squared 
modulus of the ground state excitonic wave function for 
three different dielectric environments, for a InAs tube of 
diameters = 20 nm (Fig.[S]) and D = 100 nm (Fig. [5]), 
for both hollow (top) and filled (bottom) tubes. These 
arc two relevant cases, since the former exhibits a geo- 
metrical confinement caused by the small circumference, 
whereas the latter falls fully in the 2D regime without 
any confinement. 

As shown in Fig. [5l the wave function for small tubes 
is distributed over all the circumference, best visible in 
the homogeneous case [ee = es)- Reducing the screening 
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FIG. 4: (color online) Ground state energy of the exciton for different values of the diameter D as a function of the dielectric 
contrast es/es, with ec ~ (ec ~ ^s) for hollow (filled) QTs. The shaded area indicates the limit of excitons in a 2D 
quantum well. Color and size of the symbols (filled and empty circles for the respective QTs) are proportional to the relative 
oscillator strength, normalized to the maximum value for each material (marked with an arrow). 



by diminishing affects the wavefunctions only weakly, 
leading to a slightly increased localization, both of the 
hollow (upper panels) as well as of the filled case (lower 
panels). On the other hand, both dielectric configura- 
tions lead to similar wavefunctions, reflected in the ener- 
gies reported in Fig. 21 too. 

For large tubes of Fig. [S] the wavefunction is no more 
distributed over all the circumference, but well localized. 
Therefore the curvature of the tube has no effects on the 
exciton for larger diameters, making it fully 2D. Again, 
changing the dielectric configuration, by diminishing es 
as well as by going from hollow (top panels) to filled (bot- 
tom panels) tubes, is changing the wavefunctions only 
marginally, while the respective energies are very sensi- 
tive to it (see Fig. 



Therefore, while the diameter has a definite influence 
on the dimensionality of the excitonic states, changing 
the dielectric configuration amounts to modulating the 
mean screening with nearly no effects on the wavefunc- 
tion, shifting only the energy. 



IV. CONCLUSIONS 

We have studied theoretically the excitonic properties 
of semiconductor QTs, focusing on the influence of their 
dielectric environment and its interplay with structural 
parameters. We find that, due to the strong increase of 
the electron-hole interaction and ensuing very large ex- 
citonic binding which is possible in these structures, the 
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FIG. 5: (color online) Wave function of the exciton ground state in InAs for tubes of diameter D = 20 nm and three different 
values of es- Top row (empty circle): hollow QTs, ec ~ es. Bottom row (filled circle): filled QTs, ec = es- For each panel, an 
inset on the left (bottom) shows the cut of the total wave function at y = (a; = 0). The i-axis extends from —Rtv to Rtv. 




FIG. 6: As in Fig. [5] for tubes of diameter D = 100 nm. 



spectral properties of excitonic absorption are strongly 
dependent on geometrical parameters and dielectric en- 
vironment, with energies well below the energies of the 
dielectrically homogenous case which is always in the 2D 
regime for typical parameters. Calculations have been 
performed for InAs, GaAs and InP. The low gap material 
InAs shows a peculiar behavior, since in the investigated 
systems the exciton binding energy is a substantial frac- 



tion of the gap. The very large binding energies, their 
tunability in a wide range, and the large sensitivity of 
the excitonic response to the dielectric medium, point to 
perspective applications of these systems. 
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Appendix: Derivation of the Coulomb interaction in 
a tube 



The inner radius a ~ R — 5 and the outer radius 
b = R+6 divide the space into three regions: core {p < a), 
shell {a < p < b) and environment {p > b) with dielec- 
tric constants ec, eg, es, respectively (see Fig. [T|). The 
electrostatic potential at point r induced by an electron 
localized in the shell, i.e. with a < p' < b, screened by 
es has to obey the Poisson equation in cylindrical coor- 
dinates (with charge e = 1): 

Air 

VlV^ir, r') = 5{p - p')5{(p - 4>')5{z ~ z'). (A.l) 

Here a indicates one of the three possible regions of the 
test charge: core (C), shell (S*) or environment {E). Eq. 
()A.1[) is solved by the ansatz 



27r2 



E 



dkcos{k{z - z'))gm,a{k,p,p'), (A. 2) 



where gm^aik, p, p') is the solution of the radial Poisson 
equation in each region a 



Id/ dgra.a 



pdp\ dp 



- (k^ + ^ ]gm,a = -—S{p- p') 



(A.3) 

and can be written as a linear combination of the solu- 
tions of the homogeneous Laplace equation, i.e. modified 
Bessel functions of the first kind, /^(gp), and the second 
kind, Km{qp), with the following properties)^ 

lim Irn{x) ~ 0, lim Im{x) = oo, (A. 4a) 

lim Km{x) = OO, lim Km{x) = 0. (A. 4b) 

a:— fO a;— >oo 

Imposing that limp_j.oo ,9)ri,£;(^p) = and the finiteness 
of hmp^o gjn,c{kp), wc have 



gm,c{p) = Anlm{kp) 



(A.5a) 



7<^5 (p</>) = B<J>Ukp<^>) + C<'>K„,{kp<'>) 



gm,Eip) = AnA'm(fcp), 



(A.5b) 

(A.5c) 



where 7^(5" are no Green's functions, but solutions of 
the (homogeneous) Laplace equation, from which we con- 
struct the solution gm.sip, p') of Eq. (jA.3|) in the follow- 
ing. Wc define — min[p, p'] and p^ = max[p, p']. 



Matching components of fields E\\ and D±^ at the inter- 
faces is equivalent tc4i 



< dgm,c 



dp 



> dg„i^E ^ 'ni,S 



dp 



dp 



at p = a, 

(A.6a) 

at p ~ b. 
(A.6b) 



To determine the last two unknowns we use the symmetry 
of the Green's function gm,sip,p') with respect to the 
exchange of p and p' making;^ 



gm.sip, p') = 7m,s7m,S' 

and normalization 



< , dj>s{p) 

lm.,siP) 



47r 



dp dp esp- 



(A.7) 



(A.8) 



Defining a ~ es/^c ^md l3 = es/^s, and the quantities 
K,n{kp)I,n{kp)' 



SLikp) = 



{kp)K,n{kp)' 



p=i 



.In \ K,nikp) 



TLikp) 



Im{kp) 



(A.9a) 
(A.9b) 



p=i 



Ura = T^{S^^ - a){Sl - /?) - SI^T^ia - l)(/3 - 1), 

(A.9c) 

where I{x)' = dl{x)/dx and I = a,b indicates the inner 
and outer radius of the cylindrical shell, the coefficients 
in Eq. (IA.5I) are given by 



Am 
^ m 

r> 



— ar,^(5,°„-l)(5,''„-l)/C/„ 

'"r,^(^^ 



-T-ASl~l){S:„^a)/U„ 



1' 



^SUSL~l){a-l)/U„ 
es 



Ob 

1 



1 



(A.lOa) 

(A.lOb) 

(A.lOc) 

(A.lOd) 

(A.lOc) 
(A.lOf) 



In particular, for two charges localized at the same dis- 
tance p = p' = R from the center 

gm,s{k.R,R) = [B<Im{kR) + C<Km{kR)] x 
[B>I,AkR)+C>K„,ikR)] 
An 



es 



Bl 



B' 



I,n{kR)K„,{kR), 
where we defined for clarity 



(A.ll) 
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f>< _ f£ f>< 
47r 



^ rn 

B> 



47r /™(fci?) 
ImikR) 



KmikR) 



(A.12a) 
(A.12b) 
(A.12c) 



Hence, taking Eq. (jA.2p in the special case p = p' — R, 
gives the Coulomb potential for two particles localized in 
the shell on a cylindrical surface of radius R 



For ec = = this reduces to the usual form l/es\r~ 
r'\ in cylindrical coordinates,— while for tc = (-s 7^ £b 
Eq. (|A.13|) reproduces the result of R ef . IH Note that V 
is scalable, since all arguments in Eq. (jA.13[) are products 
of lengths and momenta and thus dimensionless, only 
the measure dk of the integral is reciprocal in length. 
The latter one scales with the effective bohr length a|; — 
{es / p)0.053 nm and therefore V itself with the effective 
Hartrec Ha*= (/i/e|)27.21 cV. 



V{r,r') = — 

ires ^ — ' 



rn— — 00 



Bl 



CI 



Bl 



Cl 



l,n{kR)K„,{kR) cos(fc(z - z'))dk. (A.13) 
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